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Abstract 


This  report  examines  the  radio  propagation  model  for  narrow  and  long  tunnels. 
Modal  analysis  is  used  to  model  the  path  gain  in  2-D  and  3-D  rectangular  tunnels  and  the 
coupling  loss  of  L,  T  and  cross  tunnels.  Modal  attenuation  is  determined  by  the 


wavelength  X  and  given  tunnel  dimension  through  the  factor 


X2, 


,  where  vv  is  the 


width  of  tunnel,  by  the  mode  number  and  dielectric  constant  of  the  walls.  In  oversized 
tunnels,  many  modes  may  propagate,  however  higher  modes  have  high  attenuation 


constant  so  that  a  few  lower  modes  are  dominant  at  long  distances  from  the  transmitter. 


The  antenna  polarization,  giving  the  least  path  loss  in  3-D  tunnels,  is  found  by  comparing 
attenuation  constants  of  both  vertical  and  horizontal  polarization.  For  the  propagation  into 
branches  off  the  main  tunnel,  mode  coupling  at  the  discontinuities  is  obtained  by  hybrid 
ray-mode  conversion.  By  accounting  for  the  mode  diffraction  at  each  comer,  we  compute 
the  coupling  loss  into  cross  junctions,  the  T-junction  and  L-bends.  Because  the  cross 
tunnel  has  four  edges  and  they  generate  more  corner  diffracted  modal  fields  than  other 
structures,  and  have  the  strongest  coupling.  The  coupling  loss  is  found  to  increase  slightly 
with  tunnel  width. 
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1.  Introduction 


Radio  communication  in  mine,  rail  and  road  tunnels  has  been  of  interest  since  1920s. 
Early  studies  considered  the  use  of  leaky  cables,  which  are  widely  used  to  provide 
relatively  uniform  coverage  over  a  tunnel  [1],  However  the  use  of  leaky  cables  requires 
prior  access  to  the  tunnel  in  order  to  install  the  cable  and  related  infrastructure.  Free 
propagation  between  radio  transceivers  has  also  been  studied.  In  this  case  the  tunnel  acts 
as  lossy  waveguide.  Both  modal  [2]  and  ray  [3,4,5]  approaches  have  been  used  to  model 
propagation  in  tunnels,  as  discussed  below. 

For  small  tunnels  and  at  lower  frequencies,  the  tunnel  dimensions  are  only  a  few 
times  the  wavelength  of  the  radio  signal.  In  this  case  the  modal  approach  is  simpler  and 
more  accurate.  For  long  distance  along  a  straight  tunnel,  the  modal  approach  is 
appropriate  even  when  the  tunnel  dimensions  are  much  larger  than  the  wavelength.  The 
purpose  of  this  study  is  to  evaluate  the  coupling  of  such  modes  at  junctions  and  bends  in 
tunnels. 

For  the  modal  analysis  approach  to  propagation,  the  tunnel  is  treated  as  a  hollow 
wavegude  with  lossy  dielectric  walls.  Mahmound  and  Wait  [6]  studied  the  radiation  of 
dipoles  in  rectangular  tunnel  assuming  the  side  walls  to  be  perfectly  conducting.  Emslie 
et.  al.  [2]  found  approximate  mode  solutions  for  the  propagation  constant  and  the  field 
distribution  in  straight  tunnels  with  dielectric  walls.  Numerical  evaluation  of  the 
couplings  at  L-bends  and  T-junctions  were  carried  out  by  Sakai  and  Koshiba  [7]  using  the 
boundary-element  method.  Their  approach  gives  the  electromagnetic  field  distribution  for 
two-dimensional  tunnels  (analogue  to  to  parallel  plate  waveguide).  However,  this 
numerical  approach  does  not  yield  analytic  expressions  for  the  mode  coupling  that 
explicitly  contain  the  dependence  on  geometry,  frequency  and  polarization.  Mariage  et  al. 
[8]  investigated  field  variation  inside  or  near  a  tunnel  entrance  using  the  uniform  theory 
of  diffraction  (UTD)  to  explain  the  coupling  between  free  space  and  the  tunnel.  The  UTD 
approach  does  result  in  analytic  expressions.  In  what  follows,  we  make  use  of  the  UTD 
approach  to  describe  coupling  between  modes  at  tunnel  junctions  and  sharp  bends. 
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In  a  tunnel  whose  size  is  very  large  relative  to  the  wavelength,  another  approach  to 
propagation  modeling  is  based  on  ray  tracing.  Using  a  3-D  ray  launching  approach, 
Cichon  et.  al.  [3,4]  studied  the  radio  channel  characteristics  of  various  tunnels.  Chen  and 
Jeng  [5]  proposed  studying  propagation  in  tunnels  with  traffic  by  the  shooting  and 
bouncing  ray  method.  Ray  tracing  is  a  very  powerful  method  to  model  the  wireless 
channel  at  high  frequencies  when  the  entire  geometry  of  the  tunnel  is  taken  into  account. 
However,  for  long  relatively  narrow  tunnels  many  multiply  reflected  rays  must  be  taken 
into  account.  Moreover,  the  ray  tracing  approach  does  not  separate  the  effects  of 
individual  junctions  or  bends.  In  order  to  explain  the  propagation  characteristics  of  long 
and  narrow  tunnel  and  low  operating  frequencies,  we  make  use  of  modal  analysis.  We 
derived  the  modal  fields  and  their  excitation  by  source  starting  with  a  rigorous  integral 
formulate  for  a  straight  tunnel.  Coupling  at  junctions  and  bends  is  then  treated  using 
hybrid  ray-mode  theory  to  find  the  equivalent  sources  to  represent  the  mode  coupling  that 
results  from  diffraction  at  the  appropriate  comers.  This  approach  is  similar  to  study 
reflection  at  the  open  end  of  a  perfectly  conducting  waveguide  [9],  but  must  account  for 
the  angle  dependence  of  the  reflection  coefficient. 

We  start  by  considering  the  simplified  case  of  a  2-D  tunnel  (analogues  to  parallel 
plate  waveguide)  and  then  show  how  the  resulting  expressions  are  generalized  to  the  case 
of  actual  3-D  tunnel. 
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2.  Modal  Analysis  of  the  Radio  Propagation  in  2-D  Straight  Tunnels 

2.1  Analysis  of  TE  Propagation  in  2-D  Tunnels 


Figure  2-1.  Propagation  in  2-D  Straight  Tunnel 

Consider  propagation  in  simple  2-D  tunnel  shown  in  Figure  2-1.  The  side  walls 
have  the  dielectric  constant  ew  and  finite  conductivity  a  .  Suppose  a  line  source  of 

current  source  is  located  at  x  =  x  and  z  =  0  so  that  the  current  density  is 
J  =ayJ0S(x-x')S(z)  (2-1) 

The  field  due  to  the  line  source  excitation  can  be  found 

Vx£  =  -jcojuH  (2-2) 
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V  x  H  -  jcoeE  +  J 


(2-3) 


Substituting  expression  (2-1)  for  the  current  gives 

V  x  H  =  jcoeE  +  a  v  J0S{ x  -  x')S{z)  (2-4) 

By  taking  the  curl  of  both  side  of  (2-2)  and  using  (2-4),  we  obtain 

V  x  (V  x  E)  =  —  jo)u(V  x  // )  =  -jcoju(jcoeE  +  aYJ0S(x  -  .x')S{z))  (2-5) 

The  curl  operations  can  be  expanded  as 

Vx(Vx£)  =  V(V-£)-V2£  =  -V:£  (2-6) 

Since  V  •  E  -  0 ,  from  (2-5)  and  (2-6)  we  obtain  the  wave  equation. 

-V2E  +  k2E  =  -jcojuJ()S(.x  -  x')S(z)ci  v  (2-7) 

where  k  =  co^f/Ie  is  the  wavenumber  which  is  different  in  the  media. 

The  solution  of  the  equation  (2-7)  can  be  found  using  the  Fourier  representation 
of  the  modal  field  in  tunnel.  Since  E  is  seen  from  (2-7)  to  have  only  a  y-component,  the 
Fourier  representation  takes  the  form. 

E  (.v,  z)  =  —  f  K)e~IktdK  (2-8) 

2.71 

where  K  is  the  wavenumber  along  z.  Conversely,  E(x,k)  is  found  from  the  Fourier 
transform. 
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E{x,k)-\  E(x,z)eJKdz 

J  — OO 


(2-9) 


Taking  the  Fourier  transform  of  both  side  of  the  equation  (2-7), 

£[-v=£ 

Substituting  equation  (2-8)  for  Ev(x, z)  into  (2-10),  after  some  manipulation  we  obtain 


v  +  k 2 Ey \JKdz  =  J  [-  j(ouJQS{x  -  x')S(z)yjKdz  =  -ja)juJ0S(x  -  x')  (2- 10) 


e(x,k)  =  -jcofxJ  0S(x-  x)  (2-11) 

v  J 

To  satisfy  (2-11),  E(x,k)  must  be  continuous  at  x  =  x  ,  but  its  derivative  is 
discontinuous  above  the  plane  *  =  x  and  below  the  plane  in  the  absence  of  the  walls. 
The  derivative  must  satisfy  the  jump  condition. 


d_ 

dx 


E(x,k) 


=  o 


(2-12) 


The  electric  fields  are  given  by 

E*(x,k)=  E0(K)e~j^x~x^  for  x  >  x'  (2-13) 


and 

E  ~  (x,  k)  =  F0(/c)e  J^x~x  *  for  x  <  x ' 


where  /?  = 


-  K2  is  the  wavenumber  along  x. 


(2-14) 


5 


(2-15) 


E(k)=°*j 

2/3 


The  fields  (2-13)  and  (2-14)  may  be  viewed  as  plane  waves  radiated  by  the  line 
source.  These  waves  are  multiply  reflected  by  the  tunnel  walls.  The  total  field  E  can 
therefore  be  found  by  summing  all  of  the  multiply  reflected  plane  waves. 

From  Figure  2-1  the  transform  E(.x,k)  for  .v  >  .v'  is  therefore 


E(x,k)  =  {e-M*-*’'  +(eJflix+x')  +e-mx*x’')  £(rV;:/i/)" 


1.3.5. 


+  )  £ ( [> )" ) ||  J„ 


2  A.-  ■ 


(2-16) 


The  first  term  in  the  brackets  represents  the  wave  propagating  directly  from  the  source 
plane  .v  =  x  to  the  receiver  plane.  The  second  term  represents  waves  that  undergo  an  odd 
number  of  reflections  in  the  tunnel  walls,  initially  propagating  upward  from  the  source 
and  downward  towards  the  observation  point,  or  downward  from  the  source  and  upward 
to  the  observation  point.  The  last  term  represents  waves  that  propagate  downward  from 
the  source  and,  upward  to  observation  point,  and  upward  thus  experience  an  even  number 
of  reflection  in  the  wall.  The  power  series  may  be  summed  in  closed  form  to  give 


'g-;/?u-/»  +2cos^(.v  +  .v,)rf>~;I/t/  +  rV/4/V/?,A~t'1 1  om 
i-rVy4/i/  1 2/3 


J „ 


(2-17) 
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Eigenvalues  of  the  Propagation  Modes  in  a  2-D  Tunnel 


Using  the  Fourier  representation  (2-8),  we  can  express  the  fields  radiated  by  a 
vertical  polarized  antenna  into  a  straight  2-D  tunnel  as 


,  v  _  1  f~  f  +  2  cos  fi(x  +  x)re'J2p“  +  T  V■,4/V*Jr~y,  coju 

vi'Y’Zj_2^J-~{  1-TV'4*  2(3 


J0  \e  Jk~clK 


(2-18) 


The  modes  are  from  the  pole  singularities  of  the  (2-18)  and  are  the  solutions  of  the 
transverse  resonance  equation. 


r2  =em 


(2-19) 


For  TE  mode  (Vertical  Polarization),  the  reflection  coefficient  is 


top  top 

p7E  _  _Pir  P  _P~Ph 

(Op  +  wp  (3  +  (3lt. 

p,  p 


(2-20) 


where  (3  =  ^Jk^-K2  is  the  transverse  wavenumber  in  air  region  and  j3H,  -  -\j£n.lcQ '  -K2 

is  the  transverse  wave. .umber  in  dielectric  wall  region. 

Because  the  fields  in  the  tunnel  are  the  sum  of  inhomogeneous  plane  waves  the 
experienced  small  attenuation  along  z,  /?  and  /?u.  will  be  nearly  real  and  Flh  will  be 
close  to  -1.  In  this  case  the  solutions  for  (3  of  the  transverse  resonance  equation  (2-20) 


will  be  close  to  —  for  n=l,2,3, 
2  cl 


In  the  view  of  the  foregoing  discussion,  let 


) 3  =  --S 
2d 


(2-21) 
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Substituting  (2-21)  into  the  left  handside  of  (2-19),  the  transverse  resonance  equation 
becomes 


T2  =  e  •/4<’/5  (2-22) 

Since  Y  is  near  -1  for  the  low  attenuation  modes,  taking  the  square  root  of  (2-22)  gives 

Yte  =  -e~j2d8  «  -(1  -  j2dh)  (2-23) 


Now  consider  the  Taylor  expansion  of  Tte  about  8  =  0.  Since 
A,  =  -1)V  +  A'  ,  the  first  two  terms  in  the  expansion  are 


r7£ 


r7£(£  =  o)+<? 


2  (A -A,) 
A,(A  +  A„) 


(2-24) 


In  (2-24),  A  =  ~  and  A„  =  J(f„.  - 1  )ktt2  + 


(  nn 3" 


V2  dj 


.  Equating  (2-23)  and  (2-24),  we  may 


solve  for  S  and  obtain  the  modal  solutions. 


S.  =  ■ 


i  +  rr/r(£  =  o) 


2/3 


A  +  A, 


AA, 


jd 


A -A, 
0.(0  *0.) 


jd 


0-0. 

0.(0 +  0.) 


A„  -  A  +  A/A,  (A  +  A, ) 


(2-25) 


In  the  denominator  of  (2-25),  (3u.r/  »  1  so  we  may  further  approximate  the  solution  as 


S 


A 


"  J 


i+L' 
P) 


jd 


i+Ji+k-0 


2JK' 1 

nn 


i\ 


(2-26) 
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For  the  lower  modes  2 dk0  »  nn ,  and  hence 


71 71/ 


S„  = 


2dkn 


nA 


jdfc  i  i 


(2-27) 


Substituting  (2-27)  into  (2-22)  gives  the  poles  of  the  integrand  in  (2-18) 


n  nn  .  nA 
fa  =7T7  +  J—T 


2d  J  Ad1 


(2-28) 


The  wavenumber  along  z  for  the  n,h  mode  is  then  given  by 

Since  the  imaginary  part  fin  is  small,  Kn  may  be  approximated 


(2-29) 


.  1  1 

f ,m  ) 

1 

1 

o 

n2nX 


\k2-{- 

[2d; 


4r/37su.-l 


(2-30) 


From  (2-30)  it  is  seen  that  since  d  »  A  ,  the  attenuation  constant  is  approximately 

proportional  j  .  Assuming  that  d  »  A ,  for  the  lower  modes  having  —  «  — ,  k 
d'  2d  A 

can  be  further  approximated  as 

M2 


®  *0  -  7 


2 


(2-30A) 


where  £„.  is  now  taken  to  be  the  real  part  of  the  dielectric  constant  of  the  wall. 
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2.3  Modal  Fields  for  TE  Polarization  in  a  2-D  Tunnel 


The  modal  representation  for  the  field  in  a  2-D  tunnel  is  obtained  by  deforming 
the  contour  of  integration  in  (2-18)  so  as  to  capture  the  poles  representing  the  waveguide 
modes.  To  do  so,  it  is  convenient  to  change  the  variable  of  the  integration  (2-18)  from  K 

to  P .  Since  k  =  -yjk^  -  ~p2  ,  differentiation  of  both  sides,  gives 

dK=  ,  dP  (2-31) 


Thus  the  integration  of  (2-18)  becomes 


£,(A,:)  =  -^2„r  F(x,.\',p)ip 

Arr 


47T 


(2-32) 


where 


|  a 


+  2c0i/?(.r+ .t')rg-'vl1  +  rV'J/V"1'"'1 

”  V*V-(52 


i-r-f 


-  ;■»/*/ 


(2-33) 


The  integration  path  in  the  complex  ft  plane  lies  around  branch  cut 


K  -  Jk{)2  -  P"  ,  as  shown  in  Figure  2-2.  The  branch  points  are  at  P  =  ±A:0 ,  the  poles  are 

located  at  p  =  0  and  (3  =  —  +  j — _  ^  .  ..  for  the  TE  polarization. 

2d  4d~yjew  - 1 
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Figure  2-2  Integration  Path  in  complex  (3  plane 


By  Cauchy’s  residue  theorem,  the  contour  around  the  branch  cut  may  be  deformed  into 
the  real  [3  axis  by  summing  the  residues  of  the  poles  crossed  in  the  deformation.  If  CRl 
and  CR2  are  two  segments  of  the  contour  for  R  — >  °° ,  it  is  easily  shown  that  the 
integration  over  CR]  and  CR-,  vanish.  Thus 

J  F(P)t/(3  =  2ti;X res{F(p)}-J“oo F((3)r/(3  (2-34) 

P. 
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Since  F(/3)  is  odd  function, 


[j{p)dp  =  0 


(2-35) 


Because  the  only  contribution  in  (2-34)  comes  from  the  residues  in  plane. 
We  have 


Ey{x,z) 


2d 


( 


s 


sin((3,rv)sin((3,rv') 


e~jy[kJ- fV" 


,-jA-  iV- 

+  S  cos(P>,.y)cos(P<i.v/)  ■  2  . = 

VV  -P,.“ 

(2-36) 
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2.4  Modal  Fields  for  TM  Polarization  in  a  2-D  Tunnel 


By  the  principle  of  duality,  we  can  find  the  magnetic  modal  fields  in  2-D  tunnel 
for  excitation  by  a  magnetic  line  source.  The  solution  is  found  from  (2-36)  by  replacing 
ju0  by  £0 .  Thus 


Hv{x,z)  = 


0L)ec 

~2d 


A  e 


X  sin(P„x)sin((3#rv') 


-;V*o2-P.2: 


V*o2-P„2  n=' 


+  X  c°s (M)c°s(P/rv') 


-j-JkJ- p7- 


odd 


V*o2-P„2 


(2-37) 


However  the  wavenumvers  must  be  found  accounting  for  the  TM  reflection 
coefficient  given  by 


P„. _ L_ 

TM  _  tOepBu-  C080E„. 

,  p 

tOEpB,,,.  (0^q£)v, 


The  first  two  terms  of  the  Taylor  series  expansion  of  f,w  about  5  =  0  give 


Eu3  +  3u 


r 


TM 


=  V™(S  =  0)+S 


2ej£-0f) 
A  ,k,A  +  Aj 


(2-39) 


Making  the  same  approximations  that  lead  to  the  expression  (2-28)  for  /?„  in  the  case  of 
TE  polarization,  we  find  that 


nn  .  n  ew 

2d+J  4rf!Ve,-l 


(2-40) 
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2.5  Mode  Attenuation  in  2-D  Tunnels 


Many  propagating  modes  exist  in  tunnel  when  the  operating  frequency  is  high. 
Because  each  mode  has  different  attenuation  constant,  only  lower  modes  are  important 
for  long  tunnels.  The  attenuation  constant  for  a  TE  mode  is  defined  by 

aTE  ( dB  /m)  =  20 \og(e Im(K" } )  (2-4 1 ) 


aIL  (clB/m) 


20Im(K„)_  20 

1 

ir 

[T2  j 

In  1 0  In  10 

_2v 

Lif  1 

lyJJ 

10  >r_ 

In  10  7^  Y 


(2-42) 


where  the  parameter  y  is  defined  by 


r'W  (2‘43) 

From  equation  (2-42),  we  observe  that  the  mode  attenuation  is  dependent  on  the  factor 
- — —  and  the  square  of  mode  index  number.  That  means  that  the  lower  modes  are 

(24 

dominant  far  from  the  transmitter,  and  that  narrow  tunnels  have  higher  attenuation  than 
wide  tunnels  for  given  operating  frequency.  Similarly,  we  can  define  the  TM  mode 
attenuation  constant  as 


a 


TM 


(dB)  = 


10 


imoV^T 


(2-44) 


In  figure  2-3  we  have  plotted  the  attenuation  constant  for  the  dominant  n=l  mode 
versus  tunnel  dimension  and  wavelength  assuming  eu.  =  6  .  For  example,  a  450  MHz 

radio  wave  propagating  in  a  4  m  wide  tunnel  has  y  =0.0069,  so  that  the  attenuation  of 
dominant  TE  mode  is  1.349  dB  per  100  m.  For  the  n=2  mode,  the  attenuation  constant  is 
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Attenuation  per  100  meter  (dB) 


5.395  dB  per  100  m,  while  for  n=3,  4  the  attenuation  constant  is  12.139  dB,  21.58  dB  per 
100  m,  respectively.  Hence  only  a  few  lowest  modes  contribute  the  received  power  far 
from  the  transmitter. 


Figure  2-3  Attenuation  constant  versus  tunnel  dimension  and  frequency 

^2 

through  the  parameter - - 

(2  df 
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Attenuation/1 00m  vs  Mode  for  450MHz,  4m  Wide  Tunnel 


Mode  number 


Figure  2-4  Attenuation  constant  versus  mode  number 
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2.6  Fields  in  a  2-D  Tunnel  for  Point  Source  Excitation 


Figure  2-5  Point  Source  radiation  in  y-z  plane 

We  have  previously  derived  the  modal  expressions  of  the  fields  excited  by  a  line 
source.  We  need  to  modify  the  expressions  for  excitation  by  a  small  dipole  (point  source). 
To  do  this,  we  first  consider  the  dipole  to  be  constructed  from  a  spectrum  of  phased  line 

sources  of  strength  with  wavenumber  ky  ,  as  suggested  in  Figure  2-5.  Thus  the  y 
271 

variation  of  the  dipole  is 


I0h5(y)  =  I0h— J^e^^dky 


(2-45) 


Now  suppose  that  the  modal  fields  Ey  for  the  phased  line  source  of  unit 
amplitude  is  written  in  the  form 


Ev(x,y,z)='ZAnE„(x)e  J',ye 


(2-46) 
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n  n 

where  An  is  the  amplitude  modal  fields,  Kn  =Jkn“-ky"  and  kn  is  wavenumber  of 

mode  in  y-z  plane.  Then  from  (2-45)  and  (2-46),  the  electric  field  due  to  the  point  source 
is 


'}k'ye~jK":dkx 


(2-47) 


When  y  is  small  component  to  z,  primary  contribution  to  the  integral  (2-47)  comes  from 
the  vicinity  of  kv  =0.  In  this  case  we  can  approximate  Kn  as 


Using  (2-48),  the  integration  in  (2-47)  is  approximated  as 


Ui 


Ey(x>y>z)~^2AHEM-Jk'ye 


2  k„ 


"  ) 


elk , 


(2-48) 


(2-49) 


At  large  distance  down  a  tunnel,  the  height  y  of  the  observation  point  will  be 
about  to  the  same  as  that  of  source,  so  that  we  may  set  dipole  y  =  0  in  (2-47)  to  obtain 


Ey(x’z)=;-^SAnEn(x> 
2n  n 


-iM 


J7ITZ 

e  “k"  dk, 


(2-50) 


The  integration  of  (2-50),  can  be  evaluated  explicitly  as 


(2-51) 
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Hence  we  have  the  modal  field  from  the  point  source  at  y  =  0 , 


EA*’Z)ss-^l=LeJ *'ZAnE„(x)e  JK"Z 


s 


(2-52) 


nz 


For  the  TE  polarization,  the  modal  fields  sum  is  given  by 


JKnZ  _  _ 


2  d 


£  sin(p;rY)sin((3/rv')- - +  £  cos(P,i^)cos((3/I/)-e 


n-odd 


K. 


(2-53) 


With  the  help  of  (2-52)  and  (2-53),  the  electric  field  at  the  receiver  point  due  to  a  dipole 
is 


E  (x,  z )  =  -e 


_  Jj  Wfl0 


V2txz  2d 


£  sin(P/«JC)sin(P«-^0— - +  £  cos(P„x)cos(P„jf')- - 

=even  n=odd  ^/i 

(2-54) 


where  Kn  =  ^/k02  "Pn2  and  Pn  =  ^7  + 


nA. 


2d  \d2JeT  -1 


for  TE  mode 


For  convenience  in  evaluating  the  communication  channel,  it  is  desirable  to 
express  the  results  in  terms  of  path  loss  between  isotropic  antennas,  or  its  reciprocal 
called  the  path  gain.  In  order  to  make  this  inteipretation,  recognize  that  the  power 
radiation  by  the  Hertzian  dipole  is 


PT  =  — RrI2 

1  ?  r 


!i2 

? 


271  f 


■n 


V 


=  40^(1- h)2 
A 


(2-55) 
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where  r|  =  1  20tt  is  the  impedance  of  free  space.  The  effective  power  radiated  by  the 
dipole  in  z  direction  PERP  =  G(lPr  where  the  antenna  gain  Gd  for  a  Hertzian  dipole 


=  —  .  Thus  the  equivalent  amplitude  of  the  dipole  is 


ERP 


40^(/„-A)2=^ 

A"  Gd 


(2-56) 


Using  (2-55)  and  (2-56)  the  dipole  current  amplitude  equivalent  to  radiation  of  PERP 
watts  by  an  isotropic  source  is 


/n-/l  = 


\n  j 


P ERP 

40  G,, 


(2-57) 


If  the  effective  radiated  power  is  PERP  =  1  Watt,  then 


V'«  = 


a  5 


yTZj 


(2-58) 


The  modal  fields  for  PERP  =  1  Watt  Hertzian  dipole  is  therefore  given  by 


EY(.x,z)  =  -eJ* 


V 1 207tc 


2nd 


X  sin(p„.v)sin(pfrv') 


+  X  cos(p,rv)cos(P,rv') 

n  -odd 

(2-59) 


where  /?„ 


U7l  .  72/1 

2d+J(2<lfJe„-l 


for  Vertical  Polarization. 


The  power  received  by  a  point  source  can  be  expressed  in  terms  of  the  field  at  the 
receiver  via 


Prx  ~Grx 


(2-60) 


X2  ^|£(.x,z)|2 

4rc  2ri 

J 

For  an  isotropic  receiver,  GRX  =1,  substituting  (2-59)  into  (2-60)  gives  the  received 
power  for  PERP=l  Watt,  which  is  equivalent  to  the  Path  Gain  for  the  radio  link  (path  gain 
is  the  inverse  of  the  path  loss). 

In  Figure  2-6  and  2-7,  we  have  plotted  the  path  gain  as  a  function  of  distance 
between  isotropic  antennas  for  propagation  of  a  4  m  wide  2-D  tunnel  at  frequency 
ranging  from  30  MHz  to  1.8  GHz.  The  dielectric  constant  of  side  walls  is  taken  to  be 
£ r  =6  and  the  conductivity  o  =0.01.  The  Path  Gain  is  averaged  over  the  width  of  the 
tunnel  and  normalized  to  the  received  power  at  10  m  from  the  transmitter.  At  small 
separation,  the  path  gain  is  nearly  that  of  free  space,  which  for  isotropic  antenna  as  at  10 

(  X  f 

n:  separation  is  -  .  In  Figure  2-6,  the  antenna  polarization  is  vertical  and 

l  40tt  J 

transmitter  is  located  at  the  mouth  of  tunnel  and  2  m  from  the  tunnel  side  wall.  Because 
of  the  nairow  width  of  the  tunnel,  at  100  MHz  and  below  the  signal  decreases  rapidly 
with  distance.  The  path  gain  for  horizontal  polarization  of  the  antenna  is  shown  in  Figure 
2-7.  Because  the  attenuation  constant  of  the  TM  (Horizontal  Polarization)  modes  is 
greater  than  for  TE  (Vertical  Polarization)  modes  for  given  tunnel  width,  we  see  that  the 
signal  decreases  more  rapidly  for  horizontal  polarization  than  for  vertical  polarization. 
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Normalized  Path  Gain  of  2-D  4m  wide  Tunnel  for  Horizontal  Pol. 


Distance  (m) 


Figure  2-7  Normalized  Path  Gain  of  the  4m  wide  2-D  Tunnel  for  various  frequencies 

(8  =6,  a  -  0.01 ) 


Horizontal  Polarization 


3.  Mode  Coupling  at  L,  T,  and  Cross  Junctions 


3. 1  Mode  diffraction  analysis  of  coupling 


Figure  3-1  Propagation  in  2-D  Cross  Junction  Tunnel 


In  order  to  illustrate  the  approach  used  here  to  couple  into  side  tunnels,  consider  a 
source  located  in  one  arm  of  a  cross  junction  of  a  2-D  tunnel,  as  shown  in  Figure  3-1. 
Since  k0d  » 1 ,  the  modes  excited  by  the  source  in  the  left  arm  will  propagate  to  the  gap 

and  couple  with  little  attenuation  into  right  arm.  They  will  also  couple  into  the  upper  and 
lower  arms  as  a  result  of  the  illumination  in  the  aperture  of  the  arms.  As  a  first 
approximation,  the  illumination  in  the  aperture  is  given  by 

=cos(k;,c)+  y'sin(K„z).  For  the  lower  modes  in  the  left  arm,  K,(  is  close  to  k0, 
and  hence  is  much  larger  than  the  transverse  wavenumber  (3,,  of  those  modes.  As  seen  in 
the  cross  arms,  the  wave  number  k„  in  the  aperture  corresponds  to  the  transverse 
wavenumber  P„  in  the  cross  arm.  Thus,  the  aperture  fields  excite  only  the  very  high 
order  modes  in  the  lower  and  upper  arms  that  have  transverse  wavenumber  close  to 
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k0.  Since  these  modes  have  very  high  loss,  they  cannot  contribute  to  the  fields  at 

significant  distances  along  the  upper  and  lower  arms. 

In  view  of  the  foregoing  argument,  coupling  to  low  modes  must  occur  through  a 
different  mechanism.  Felsen,  in  studying  reflection  of  modes  incident  on  the  open  end  of 
a  waveguide,  suggested  such  a  mechanism  [9].  Viewing  the  incident  modes  as  composed 
of  multiply  reflected  plane  waves,  the  coupling  to  the  reflected  modes  resulted  from 
diffraction  at  the  edges  of  the  open  end  waveguide.  Each  plane  wave  excites  diffracted 
rays  that  are  multiply  reflected  in  the  waveguide  walls.  By  summing  the  diffracted  ray 
contributions  and  using  the  Poisson  sum  formula,  it  is  possible  to  translate  the  diffracted 
rays  into  a  modal  sum,  with  the  amplitude  of  the  modes  given  by  the  diffraction 
coefficient  evaluated  at  the  mode  angle.  The  diffracting  corners  and  the  rays  that  excite 
the  cross  tunnel  are  shown  in  Figure  3-1.  This  approach  is  especially  will  suited  to  the 
case  when  k0d  » 1 ,  as  we  have  found  necessary  for  propagation  over  a  significant 

distance  in  a  tunnel.  The  multiply  reflected  rays  originating  at  the  upper  right  hand  corner 
of  the  upper  arm  in  Figure  3-1  are  shown  Figure  3-2.  In  this  figure  we  have  reoriented  the 
coordinate  system  so  that  z  lies  along  the  axis  of  the  cross  tunnel. 

The  field  received  at  a  distance  z  along  the  tunnel  is  the  sum  of  the  fields  earned 
by  the  multiply  reflected  ray  produced  by  diffraction  at  the  comers  of  the  tunnel.  The 
multiply  reflected  rays  arriving  at  the  receiver  may  be  viewed  as  coming  from  multiple 
images  of  the  edges  in  the  tunnel  walls.  The  image  array  for  a  single  line  source  at  the 
right  hand  edge  is  shown  in  Figure  3-2.  Similarly,  the  left  hand  edge  gives  the  diffracted 
ray  contributions  that  appear  to  come  from  the  images  shown  in  Figure  3-3. 

Each  mode  incident  on  the  junction  will  generate  a  set  of  sources  of  the  type 
shown  in  Figure  3-2  and  Figure  3-3.  The  source  corresponding  to  the  right  hand  edge  will 
have  strength  »0  that  >s  determined  by  the  phase  and  amplitude  of  the  incident  modal 

field,  and  a  pattern  function  D(0;i,0,„ ).  Here  0,,,  is  the  angle  of  incident  plane  wave 
comprising  the  mode  incident  on  the  junction,  and  0,,  is  the  angle  of  the  diffracted  ray. 
Similarly,  the  left  hand  edge  will  have  amplitude  m,  and  pattern  function  D(Qn,Qin ).  The 
explicit  expression  for  Uq  ,  and  D(0,i,0()l)  will  be  discussed  later.  When  uQ ,  u (  and 
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D(0/(,0,;1)  are  known,  then  forTE  polarization  the  field  Ex  at  the  receiver  location  (x,z) 
is  given  by 


EyU,z)=  E  E  »,r(e„)'’£>(e„,0j. 


1=0,1  «=-oo 


(3-1) 


rn  =  ^7  2  +  (.v  —  (4/?  ±  l)c/)2 


(3-2) 


0„  = tan  1 


x-(4n±l)cl 


(3-3) 


where  the  +  sign  is  applies  the  right  hand  edge  (i=0),  the  -  sign  is  for  the  left  hand  edge 
(i=l).  In  (3-1)  ~  (3-3),  n=0  refers  to  the  direct  ray  and  n  =  ±1,±2,±3,  ...  refers  to  rays 
coming  from  the  images,  as  indicated  in  Figure  3-2  and  3-3.  The  index  p  >  0  gives  the 
number  of  reflections  undergone  by  the  ray  corresponding  to  the  particular  image.  The 
values  of  p  are  listed  in  Figure  3-2  and  3-3,  and  expression  for  p  in  terms  of  n  will  be 
given  in  the  subsequent  discussion. 

In  the  case  of  a  waveguide  with  metallic  walls,  |T|  =  1  so  that  only  the  sign  of 

(r)/;  needs  to  be  accounted  for.  In  the  case  of  dielectric  waveguide  walls,  (r)p  is  a 
function  of  the  number  of  reflections  through  the  glancing  angle  0„ ,  and  through  the 
exponent  p.  In  order  to  generalizing  Felsen’s  method  to  the  case  of  dielectric  walls,  it  is 
necessary  to  consider  [p(0„ )]/;  as  a  function  of  ray  direction  so  that  the  pattern  function 

of  source  becomes  [r(0„  )p  D(0/; ,  0(„ ) .  In  this  way  the  field  from  each  image  accounts 
for  the  multiple  reflection. 
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3-2  Ray-Mode  conversion 


For  propagation  over  a  significant  distance  along  the  tunnel,  many  reflections 
must  be  accounted  for  in  the  sum  (3-1).  A  more  convenient  representation  of  the  field  in 
this  case  is  found  by  converting  the  ray  sum  into  mode  sum,  as  discussed  below  in  terms 
of  Poisson  sum  formula  [10]. 

A  periodic  function  f(x)  with  period  4d  can  be  represented  in  form 


7W=  x/(-*+4^) 

fl  —  OO 


(3-4) 


where  f(x)  is  some  function  defined  on  the  infinite  interval.  As  applied  to  the 
waveguide  problem  f(x)  is  the  y  component  of  electric  field  inside  the  waveguide 
-d  <x<d  ,  and  is  the  analytic  continuation  of  the  field  for  the  remainder  of  the  period 
4d  ( —  2 d  <  x  <  - d  and  d  <  x  <  2d ).  The  Poisson  sum  formula  expresses  this  function  as 
the  sum 


.nm 

/W ^  ” 

c  m=~ <*> 


where  Cm  is  the  Fourier  transform  of  f{x)  given  by 


(3-5) 


.nm 

Cm  =IL/M*  Jld  dx  (3-6) 

If  the  reflection  coefficient  T  =  1 ,  then  we  can  interpret  f(x)  as  being  the  y 
component  of  the  electric  field  due  to  a  line  source  inside  the  waveguide,  in  which  case 
the  term  f(x  +  And)  represents  the  field  due  to  one  of  the  images  of  the  line  source  in  the 
waveguide  walls,  as  shown  in  Figure  3-2  for  a  line  source  at  d,  or  in  Figure  3-3  for  a  line 
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source  at  -d.  When  the  waveguide  walls  are  not  perfectly  reflecting,  the  definition  of 
/(.v  +  And )  is  modified  to  include  reflection  coefficient,  as  discussed  below. 

By  grouping  terms  in  (3-5)  for  m  =  ±|//;| ,  the  sum  (3-5)  is  expressed  in  terms  of 

the  waveguide  modes.  Note  that  the  m  =  0  mode  in  the  dielectric  waveguide  would  be 
highly  attenuated  and  is  ignored  here.  According  to  Poisson  sum  formula,  the  total  field 
Ey(x,z)  in  the  waveguide  is  given  by  the  modal  expression 


£v(.v,z): 


1 

2d 


~  r  +  C 
^  ^  m  _  2Zz!!L 


COS 


nm 

2 ~d 


-x 


m=l 


C  -C 

^  n ;  —h 


■sin 


(  nm 

Id 


JJ 


(3-7) 


The  coefficient  Cm  is  given  by  the  Fourier  transform  of  an  equivalent  line  source  field 

/(.v). 


3.2.1  Evaluation  of  Cm  for  diffraction  at  the  left  hand  comer 

The  choice  of  the  equivalent  line  source  field  can  be  understood  in  terms  of 
Figure  3-3,  where  we  have  shown  the  original  line  source  at  -d  and  its  images  at  +3d, 
+7d,  ...  and  -5d,  -9d,  ....  The  images  on  the  right  are  indexed  by  n  =  1,  2,  3,  ...  and 
those  on  the  left  by  n  =  -1,  -2,  -3,  ...  .  Each  image  represents  the  ray  that  is  reflected  p 
times  in  the  waveguide  walls,  where  p  =  2n-l  for  n>0  and  p  =  -2n  for  n<0.  The  individual 
terms  in  the  ray  sum  (3-4),  which  represent  the  multiply  reflected  rays,  will  have  the 
values 


f(x  +  4nd)  =  zq 


[r(0„)]'’£>(e„.e,„) 


(3-8) 


where 
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(3-9) 


rn  =  -^/z2  +  {(.*  +  d)-4;ic/}2 
and 

0;,  =  tan  _1 

Here  ii\  is  the  amplitude  of  the  field  incident  on  the  edge  and  the  diffraction  coefficient 
D(0„  ,0,„ )  gives  the  pattern  function  of  the  line  source. 

By  replacing  x+4nd  in  (3-8)  ~  (3-10)  by  x,  it  is  seen  that 

/M = «i[r(e)]p(>)  D(e„ ,  e,„  )e-^~  <3-1  d 

Vr 


(jc  +  d)  -  4nd\ 


(3-10) 


where 


r  =  i/z2  +(x  +  r/)? 


(3-12) 


and 


0  =  tan 


r  ^ + ci  ^ 


(3-13) 


In  (3-11)  the  exponent  of  the  reflection  coefficient  is  shown  by  the  step  curve  in 
Figure  3-4 

In  order  to  carry  out  the  integration  (3-6),  we  write  Tp^  in  the  form 
[r(0)]pW  =e/7(v)lnr(0)  (3-14) 


31 


If  we  now  express  T  as  a  function  of  0  and  assume  that  z  is  large  so  that  sin  0  «  1  for 
the  lower  values  of  \n\,  then  forTE  polarization  we  have 


sin  0-yeM,  —  cos “  0 
sin  0  +  a / 8 -cos2  0 


2sin  0 


(3-15) 


We  further  approximate  In  f(0) 


In  r(0)  =  7’7l  +  In  1 


2sin  0 


2  sin  0 


£«•" 1 


(3-16) 


Using  the  approximation  in  (3-14)  gives 


[r(0)]/j(v) 


p(x)  jn — —  (  i)/?(~v)  exp  2smBp(.r)  =  -sgn(.v-2</)exp  - 

^  V£vr-lJJ  LVen-l  J  [V 

(3-17) 


Figure  3-4  Variation  of  the  reflection  coefficient  index  p(x)  for  the  left  hand  edge 


Note  that  (3-17)  is  a  generalization  of  the  approach  used  by  Felsen  [9]  for  metallic 
waveguides  where  F  =  -l  so  that  rp^  - -sgn(x- 2d).  Using  (3-1 1)  and  (3-14)  it  is 
seen 


C  = 

^  m 


"Co 11 1  D(0’  )sSn(*  -  2d)  exp 


-  2  sin  0 


p{x) 


,-jkny[7+{x+df 


jz2  +  (x  +  d)~ 


.mn 

e~JuXdx 


(3-18) 


The  function 


exp 


-  2  sin  0 


in  (3-18)  has  discontinuities  due  to  the 


definition  of  p(x)  at  x  =  ±  2d,  ±  4d,  ±  6d,  ....  However,  since  |9|  is  small  for  -  »  d ,  at 

the  discontinuities  closest  to  d,  the  jump  in  the  exponent  will  be  small.  For  discontinuities 
further  from  d,  the  exponential  is  already  small  so  that  it  will  give  small  contribution  to 
the  integration.  We  therefore  approximate  p(x)  in  (3-18)  by  the  ramp  function 


P{*)  = 


x  +  d 

—  .  This  approximation  is  shown  by  the  das/ied  lines  in  the  Figure  3-4.  The 


choice  of  this  ramp  function  is  made  to  in  order  to  simplify  the  asymptotic  evaluation  of 

-  2sin  0 


the  integration  for  Cm  .  The  variation  of  the  term  exp 


^  1 


PM 


in  (3-18)  using  the 


exact  and  approximate  forms  of  p(x)  are  plotted  in  Figure  3-5.  In  (3-18) 
(x  +  d) 


sin  9 


V- 2  +(-v  +  d)2 

Using  the  foregoing  expressions  and  the  approximation  form  for  p(x) ,  we  define 
the  functions 


nm 


(x  +  d)2 


q(x)  =  k0Jz?+(x^^2  +-—X-J-  _  _ _ 

d  z~  +  (x  +  d)z 


(3-19) 
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Figure  3-5  Linear  approximation  of  reflection  constant  ex 


P 

v 


PW 


for  the  diffraction  of  left  hand  edge,  where  0  is  given  by  (3-13),  z  =  100  m  and  d  =  2  m 


(3-20) 


Then 


C„ 


it [  sgn(.v  -  2 d)g {x)e  Jq'x  V/.v 


(3-21) 


To  simplify  the  integration  we  change  variable  to  T|  =  .v  +  cl ,  in  which  case 


qiy\)=k^z2  +ri2  +^-{r]-d)-j 


dJew  -lVz2  +rl2 


sM  = 


D(e,e, 


(3-22) 


(3-23) 


Cm  =  sgn(ri-3rf),?(r|)e  Jq{r])dr\ 


(3-24) 


3.2.1. a  Evaluation  for  m<0 


For  m  <  0  and  for  z  »  >  the  saddle  point  of  g(r|)  will  lie  in  the 


region 


T]  >  3 d  .  Then  the  integration  can  then  be  written  in  a  form  that  allows  the  separate 
evaluation  of  the  saddle  point  and  end  point  contributions.  This  form  is 


Cm  =  ~u\  t  sM®  Jq^dx  -  2tt  sOl)*  Jq{r])dx\ 


(3-25) 


For  the  region  of  interest  z»  cl  »X,  the  saddle  point  r|5  is  obtained  from  the 
solution  of  q'(r\)  =  0  .  The  derivation  of  c/(r|)  can  be  approximated  by 


//  \  M  mn  ■  2r\ 

g(Ti)«  — 7  = 

z  2d  dzJew- 1 


(3-26) 


and  the  saddle  point  is 
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It  is  seen  from  (3-27)  that  for  m<0  and  z  »  >  the  value  of  Rer|s  will  be  greater  than 

3d. 

Consistent  with  the  foregoing  approximation  we  find  that  <r/(r|s ),  g (r| y )  and 
g"(r\s )  are  given  by 


/  \  .  n  T  mn  mn  .  n  " 

2d  2  zd^jew- 1 


/H7t 

~2d 


( 


mX  .  mn 

—  +  J  — 

4  cl 


M 


d-k0-Jeu.  - 1 


7J 


mn 


*(»!,)  = 


D(0,„9,„)  D(9,„.e,„) 


+t 


VI 


where 


0m  =  tan" 


=  tan 


4r/ 

V  7 


and 


(3-28) 


(3-29) 


(3-30) 


(3-31) 


The  saddle  point  integration  Li  1]  for  Cin  with  m<0  is  obtained  as 


JI*(n)*  J(lin)dr\ 


2n 


nm 


=  if«(e,„.e„,v"yr‘M 

V  Jko 


mX  .  niTt 
4d  J  d2k02^tix~  1 


./?ITC 


<?  4-\/AD(0m,0/;i)«?  JK"'zeJ  2 


.  J- 


.mn 


(3-32) 


Note  that  the  coefficient  of  z  in  the  exponent  is  the  same  as  the  value  of  Km  given  by  (2- 
30A)  for  the  mode  solution  in  a  straight  tunnel  under  the  condition  that  d  »  X . 

The  end  point  contribution  to  Cm  is  written  in  the  form 


lend  M{tr\ 


(3-33) 


2  / 

If  z  »  d  A  ,  then  the  saddle  point  will  be  far  from  the  end  point  and  we  can  evaluate  (3- 

33)  asymptotically  by  approximating  the  exponent  of  integrand  in  (3-33)  near  the 
endpoint  as 


q(r\)  =  q(3d )+  (t|  -  3d )q{3d ) 


(3-34) 


With  the  help  of  (3-34),  the  end  point  integration  becomes 


len,  «  giidy^f^e-^^dn 


(3-35) 


Letting  x  =  j(r\-3d)q'(vl),  we  have 


sM  e-M  f°e-Tdr  =  sM  e-jc,M 

jq'{3d)  J“  jq'{3d) 


(3-36) 
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where 


cj( 3d)~  k0z  +  nm 
nm 


2d 


D 


sM=- 


-l  3r/  Q 
tan  — ,  0 ... 


D{  O,0,J 


("~  +  9d  2 ) 


(3-37) 

(3-38) 

(3-39) 


Thus  the  end  point  integration  in  Cm  is  obtained  as 


g{n)e~jciir])dr]  »  2rfD(°»Qf>.)g-y(V+>»m) 
jnmd  z 


(3-40) 


With  the  results  of  (3-32)  and  (3-40),  C...  for  m  <0  becomes 


cm  “  ~“i 


.71 


.  fUTl 


i1*4XD(am,siy  ^  4t/D(°-°J,-A,- 

jmn 


(3-41) 


For  z  »  >  then  ^  and  the  second  term  in  (3-41)  coming  from  the  end 

point  integration  is  small  compared  to  the  first  term  due  to  the  saddle  point.  We  therefore 
omit  the  end  point  term  from  later  analysis. 
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3.2. 1  .b  Evaluation  for  m>0 


For  m  >  0  the  saddle  point  lies  in  the  region  r\  <  3d  ,  and  we  may  write 

Cm  = 11 1  {d,  8  Ol  )e  ~jci^dx  -  2j^  j  (3-42) 

The  evaluation  of  saddle  point  contribution  is  like  that  for  m<0,  and  leads  to 
expression  (3-32).  Using  approximation  employed  to  evaluated  (3-33),  the  end  point 
contribution  in  (3-42)  is  the  negative  of  (3-40),  or 


5d  -  jmny/z 

Thus  the  integration  Cm  for  m>0  is  given  by 


Cm  ~U\\ 


-J 


.71 


■*VAD(e„„e,„)e 


jun 


4dD(o,ein) 

jnm4z 


-A' 


(3-43) 


(3-44) 


The  second  term  in  (3-44),  which  cam''  from  the  end  point  integration,  is  small  compared 
to  the  first  term  provided  z  »  ■  As  in  the  case  of  m<0,  we  neglect  this  term  in  later 

analysis. 
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3.2. l.c 


Total  field  due  to  diffraction  at  the  left  hand  edge 


From  the  result  of  (3-41)  and  (3-44),  and  neglecting  the  second  term,  the 
diffracted  modal  electric  field  for  z  »  due  to  the  itl  incident  field  to  the  left  edge 
in  the  tunnel  is  obtained  as 


Ev{x,z)=ul 


2  d 


Wl-I 


X(-i)  2  D(e„,,eJcos| 

m-odd 


1  mn  ^ 

2 7' 


m 

-  E(-i)i^(e„(,e(),)sin 


( mn  ' 

UXj 

J 

(3-45) 


where  K,);  is  given  by  (2-30A).  Expression  (3-45)  is  in  the  form  of  a  modal  sum,  as  in 

(2-26),  except  that  the  transverse  wavenumber  P/)(  is  replaced  by  j  ■  However,  since 

the  imaginary  part  of  the  transverse  wavenumber  in  (2-26)  is  very  small  for  cl  »X  ,  the 
transverse  variation  in  (3-45)  is  nearly  that  of  the  modal  field. 
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3.2.2  Evaluation  of  Cm  for  diffraction  at  the  right  hand  comer 


For  the  diffracted  field  due  to  the  right  hand  edge  shown  in  Figure  3-2,  the  field 
expression  in  tunnel  is  expressed  as 


/W="o[r(e)]',(')D(e„,eJ. 


,~Jkor 

7T 


(3-46) 


Here  Uq  is  the  amplitude  of  the  incident  field  at  the  right  hand  edge, 


=  -^/z2  +  („t-r/)2 


(3-47) 


and 


,/ix 


0  =  tan  ’ 


-4 


V  -  J 


(3-48) 


The  exponent  of  reflection  coefficient  p(x)  in  (3-46)  is  shown  in  Figure  3-6. 
Using  the  same  approximation  (3-!  7)  for  T ,  the  reflection  coefficient  in  (3-46)  becomes 


[r(0)4'v)  =  sgn(x  +  2<7)exp 


-  2sin  0 

Ve»'  _1 


p(x) 


(3-49) 


With  (3-46)  and  (3-49),  C,„  is  given  by 


Cm  =  lZollo  sgn(A'+  2d)D(B,Qin  )expi 


-  2sin  0 
a/8.."1 


P(x) 


,-jk0yjz:+(.x-df 

—  e  2d  clx  (3-50) 
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Figure  3-6  Variation  of  the  reflection  coefficient  index  p(.x)  for  right  hand  edge 


As  before,  we  approximate  p(.x)  in  the  exponent  of  (3-49)  by  the  ramp  function 
,  %  \x  -  d\ 

p{x)  = J - -  shown  in  Figure  3-6  to  simplify  the  asymptotic  integration.  The  integration 

2d 

for  Cm  is  evaluated  by  same  approaches  as  we  showed  for  the  diffraction  of  the  left  hand 
edge.  Thus  Cm  for  m>0  is  obtained  as 


Cm  =  “« 


0 


.71  .mu 


J* 


4</OlO.O,„  )  _ 


jkp 


jmn 


■Tz 


(3-51) 


And  Cm  for  m<0  is  given  by 


Cm  ““O'! 


.n 

~J— 


.niK 
~ J — 


4  VXD(em,0/Jc  2  e~jK" 


-(-1 Y 


4r/D(O,0,,)  , 


jnm 


V7 


(j 


(3-52) 
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As  before,  for  z  »  we  may  neglect  the  second  term  in  (3-51)  and  (3-52). 

From  the  result  of  (3-51)  and  (3-52),  the  diffracted  modal  electric  field  due  to  the 
diffraction  at  the  right  hand  edge  in  the  tunnel  is  expressed  as 


Ev(x,z)  =  -uQe 


m- 1  f 

0 ,  i  £(-i)^£>(ew,e,Jcos| 

la  \m=odd 


mn 
2 d 


m 

-  2(-0To(e,„.e,.„)sin 

m=even 


mn  ~) 

- A' 

e~K,z 

2^  J 

J 

~J  K,„ 


J 


(3-53) 
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3.3 


Mode  Diffraction  in  2-D  Cross  Tunnels 


I 


Figure  3-7  Mode  Diffraction  in  2-D  Cross-Tunnels 


For  a  2-D  cross  tunnels  we  suppose  that  the  n'1'  mode  is  incident  from  the  left  in 
the  primary  waveguide.  The  modal  fields  excited  in  the  cross  arms  are  then  generated  by 
diffraction  at  the  four  corner  edges  (C 1  ,C2,C3  and  C4),  as  indicated  in  Figure  3-7  for  the 
fields  in  the  lower  arm.  Let  Ev(x.z)(  (  C2  =  ^dif\  +  ^dif2  be  the  field  in  the  lower  cross 

arm  due  to  diffraction  at  the  left  hand  comers  Cl  and  C2.  These  contribution  to  the  total 
field  were  derived  in  the  previous  section.  Assuming  both  the  main  tunnel  and  cross  arm 
to  have  width  2d,  this  contribution  to  the  field  is  given  by  the  modal  sum 


44 


E ?{X'Z)c\,C2  - 


j~~  \  m~ *  m  1 

» „ie  4—  E(-0~ 

£Cl  \^m=odd  m-even  ) 

+  U„,e-'2*-«.JJ  f-\  Z(-l)T!D((|),:,2,t,2)cosO,„.v>-^--  -  2(- Of  D(0,„2 .<t»„2 )sin(p„,.v>-^~ 

[w=f;rtW  m-even 

(3-54) 


Here  is  the  amplitude  of  the  field  incident  on  comer  Cl,  while  the  field  incident  on 

the  comer  C2  differs  by  the  phase  factor  e~ildK" . 

In  addition  to  the  diffraction  at  the  right  hand  comers,  diffraction  at  the  upper 
comers  generate  £v(j.z)g  C4  =  Edjj 3  +  Edij 4,  which  is  given  by 


The  term  un2  is  the  amplitude  of  the  field  incident  on  comer  C3,  while  the  field 
incident  on  the  comer  C4  is  ull2e~^2dK" .  The  various  diffraction  angle  in  (3-54),  (3-55) 

are  as  shown  in  Figure  3-7.  Also,  the  fields  unl  and  un2  represent  Evn  of  the  n,h  mode 

at  the  lower  and  upper  walls  of  the  main  tunnel.  For  line  source  excitation,  they  are  given 
by 
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£>."(J:  =  ±rf,z  =  L)  =  -^-y0i 


2  cl 


sin(P„  (±  c/))sin(P,rv/) 


-j*J' 


K, 


(3-56) 


cos((3„  (±f/))cos((3„/) 


-j*nL 


K, 


n-odd 


where  L  is  the  distance  from  the  transmitter  to  the  corner.  The  total  modal  field  in  the 
cross  arm  is  then  given  by  a  summation  over  the  incident  modes  in  the  main  tunnel,  or 


Ey  ^  (*’  z)~  X  iEy  (A’>  z)ci.C2  +  Ey  (-v>  -)c3.C4 ) 


(3-57) 


To  find  the  field  in  2-D  cross  tunnels  due  to  point  source  in  the  main  tunnel 
requires  the  inclusion  of  the  field  spreading  of  the  plane  of  the  cross  section.  For  the 
lower  modes  having  K„  close  to  the  spreading  is  account  for  by  multiplying  (3-56) 


.71 

by  the  factor  e  4 


yfko 


sj2n(z.  +  L) 


.  In  order  to  study  the  effect  of  the  diffraction  into  the  cross 


arm  on  the  received  signal,  we  have  compared  the  received  power  as  a  function  of 
distance  along  a  straight  tunnel  to  the  received  power  as  a  function  of  distance  along  a 
cross  arm  located  L  =  500  m  from  the  source.  Because  of  the  various  phase  terms  in  the 
sums,  the  fields  show  considerable  local  variation  over  a  wavelength  scale.  Some  of  this 
variation  can  be  removed  by  averaging  the  received  power  over  the  cross  section  of  the 
tunnel.  Figure  3-8  shows  the  average  received  power  as  a  function  of  distance  for 
vertically  polarized  antennas  operating  at  900  MHz.  Both  tunnels  are  4  m  wide  and  the 
walls  have  dielectric  constant  er  =  6  and  conductance  a  =  0.01 .  The  1  watt  transmitter  is 
located  in  middle  of  the  tunnel  (2  m  from  the  tunnel  wall).  The  UTD  diffraction 
coefficients  proposed  by  Luebbers  [12]  have  been  used  for  the  calculation.  We  found  that 
the  coupling  loss  of  cross-junction  tunnel  is  about  37  dB  for  this  case.  Plots  of  power 
received  in  the  cross  tunnel  for  300MHz  and  100  MHz  sources  arc  shown  in  Figure  3-9 
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Receiusd  Power  (dB) 


and  3-10.  For  these  computations,  the  comers  are  located  at  L=50  m  from  the  transmitter 
and  coupling  loss  is  found  to  be  28  and  23  dB,  respectively. 


Figure  3-8  Received  Power  in  4  m  wide  2-D  Straight  and  Cross-junction  Tunnels 
900  MHz,  Vertical  Polarization,  er  =  6,  o  =  0.0 1 
The  comer  of  cross-junction  Tunnel  is  located  at  500  m. 
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Received  Power  (dB) 


Received  Power  in  4m  wide  Tunnel  for  300MHz,  Vertical  Polarization 


Figure  3-9  Received  Power  in  4  m  wide  2-D  Straight  and  Cross-junction  Tunnels 
300MHz,  Vertical  Polarization,  £,.  =  6,0  =  0.01 
The  corner  of  cross-junction  Tunnel  is  located  at  50  m. 


Received  Power  in  4m  wide  Tunnel  for  100MHz,  Vertical  Polarization 


Figure  3-10  Received  Power  in  4  m  wide  2-D  Straight  and  Cross-junction  Tunnels 
100  MHz,  Vertical  Polarization,  ew,  -6, a  =  0.01 
The  comer  of  cross-junction  Tunnel  is  located  at  50  m. 
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3.4 


Mode  Diffraction  in  2-D  T-Tunnels 


Figure  3-1 1  Mode  Diffraction  in  2-D  T-junction  Tunnels 

When  a  modal  field  is  incident  on  a  T  junction,  as  shown  in  Figure  3-11  it 
illuminates  the  two  corners  directly,  and  indirectly  via  reflection  at  the  further  wall.  For 
cl  »  A  the  reflected  field  will  be  essentially  the  same  as  that  of  two  incident  mode, 
except  for  a  reflection  coefficient  T.  Because  the  mode  angle  <J>„  are  small,  T  will  be 

close  to  the  value  for  a  plane  wave  at  normal  incidence.  Let  1  (.v,  z)  represent  the 

modal  field  in  the  lower  arm  of  the  junction  due  to  the  n'1'  mode  incident  directly  from 
the  horizontal  tunnel  in  Figure  3-1 1  into  the  upper  and  lower  comers.  This  sum  is  given 
by 
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Eydifl(x,z)  = 

Un\Z*  ~7"[  Si"1)”  D^,n\Anl)COS^mxyjK",Z 

lcl  [m=odd 

-  ^{-l)^D{^mi,(\>nl)sm{^mx)e~jK'"z  1 

m-even  J 

-  u,ae~‘  £  (  S  (- if?  D(*„i2 , 4*112  IcoslMK*-'"2"’ 

[/H=0rfrf 

-  S(-  of  0(<t>,„2 -4>„2  )sin(p„,-vV-^-<=*2rf) } 

m—even  I 


Similarly,  let  e/^2{x,z)  represent  the  modal  field  resulting  from  diffraction  at 

the  comers  resulting  from  the  incident  mode  field  that  is  reflected  from  the  wall.  This 
sum  is  given  by 

M= 

u^re-^-A^l  2(-l)T1Z)(4.„„r,4',/)cos(|i„^-'K-; 

Zcl  [m-odd 

-  S(-l)fo(4,„/,4,,/)sin(|3„.rV-jK-4 

m=even  J  (3-59) 

-U„2r S (- 1)^ o(<t>,„2r ,<t>„2" )cos(p„,.v>-^-<=+2d» 

M  [m-odd 

-  4«2r.^2r)sin(P«^-^(=+2rf)} 

m-even  \ 


In  (3-59),  T  is  taken  to  be  the  plane  wave  reflection  coefficient  for  normal  incidence. 
The  incident  fields  unl  and  un2  are  given  by  (3-56). 

The  total  field  in  the  lower  arm  is  then  given  by 


E  d,f  (x,z)  =  E  k*^1  (*.  z)+  Eydif2  (jc,  z) 


H=1 


(3-60) 
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Computation  have  been  made  for  the  received  power  averaged  over  the  cross 
section  in  the  main  tunnel  and  in  the  side  tunnels  for  a  T  junction.  The  tunnels  are  4  m 
wide  and  the  walls  have  dielectric  constant  zr  =  6  and  conductance  a  =  0.01 .  Figure  3- 
12  shows  received  power  in  straight  tunnel,  and  in  the  arm  of  a  T  junction  for  a  vertically 
polarized,  900MFlz  transmitter  located  in  the  middle  of  the  main  tunnel  (2m  from  the 
each  wall)  and  500m  from  the  T  junction.  The  coupling  loss  into  the  T  junction  tunnel  is 
about  42  dB  for  this  case.  The  received  power  for  300  MHz  and  100  MHz  sources  are 
shown  in  Figure  3-13  and  3-14  when  the  sources  are  located  at  50m  from  T  junction.  The 
coupling  loss  for  theses  two  frequencies  is  about  31  and  29  dB,  respectively. 
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Received  Power  in  4m  wide  Tunnel  for900MHz,  Vertical  Polarization 


Figure  3-12  Received  Power  in  4  m  wide  2-D  Straight  and  T-Junction  Tunnels 
900MHz, Vertical  Polarization  ,  £„,  =  6, a  =  0.01 


The  comer  of  T-Junction  Tunnel  is  located  at  500  m. 


Received  Power  (dB) 


Received  Power  in  4m  wide  Tunnel  for  300MHz,  Vertical  Polarization 
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Figure  3-13  Received  Power  in  4  m  wide  2-D  Straight  and  T-Junction  Tunnel 
300MHz,  Vertical  Polarization  ,eu.  =  6, a  =  0.01 
The  comer  of  T-Junction  Tunnel  is  located  at  50  m. 


Received  Power  in  4m  wide  Tunnel  for  100MHz,  Vertical  Polarization 
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Figure  3-14  Received  Power  in  4  m  wide  2-D  Straight  and  T-Junction  Tunnels 
100MHz,  Vertical  Polarization  ,eu,  =6, a  =  0.01 
The  comer  of  T-Junction  Tunnel  is  located  at  50  m. 
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3.5 


Mode  Propagation  in  2-D  L-Tunnels 


Ediff 


p 


Figure  3-15  Mode  Diffraction  in  L-bend  Tunnels 


When  turning  the  comer  in  an  L  shaped  tunnel,  the  diffraction  takes  place  at  the 
single  corner  shown  in  Figure  3-15.  However,  reflection  in  the  opposite  wall  produces 
additional  contributions  to  the  modal  fields.  One  contribution  to  the  modal  field  in  the 

arm  of  the  tunnel  is  directly  diffracted  field  Evlli^(.\,z)  that  is  diffracted  through  the 

angle  <\>m ,  and  caused  by  the  field  incident  at  the  mode  angle  <J>„ .  A  second  field 

Evdi^2(x,z)  results  from  the  reflection  of  the  incident  mode  at  the  wall,  seen  on  the  right 

in  Figure  3-15.  The  incident  angle  for  this  field  is  <|)„r  =  <J)„  ,  while  the  angle  of  the 

diffracted  ray  is  <)),„'  in  Figure  3-12.  The  third  field  Evcl‘^( .x,z)  from  the  reflection  of  a 
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diffracted  ray  at  the  top  wall  in  Figure  3-15.  In  this  case  the  diffraction  angle  is  . 

This  third  contribution  is  generated  by  both  the  incident  modal  field,  and  by  the  modal 
field  reflected  in  the  right  hand  wall  in  Figure  3-15.  However,  because  of  the  reduction  in 
strength  due  to  the  reflection  coefficient  for  near  normal  incidence,  we  include  only  the 
contribution  from  the  incident  normal  field.  The  modal  field  in  the  tunnel  arm  is 


E,dif  {x, z )  =  £ k^'Cv, z)  +  z)+  Eydif3 (x,  z\ 


n=\ 


(3-61) 


The  modal  diffracted  field  from  the  single  comer  Evcli^(x,z)  is  expressed  as 


E/'f'(x,z)  = 

j-  Jx  f  '"-1 

«,1  e  4ttt  EH)  2  D{$m, 0fl)cos(PlfIx)e  JKz 

[m=odd 


m 

~  EH^  )sin(P„rv>'”-/K"'; 

m=even 

(3-62) 


where  the  incident  modal  fields  un  at  the  comer  is  given  by  (3-56). 

The  fields  Evcli^2(x,z),  Evdi^3(x,z)  are  of  the  same  form  as  to  (3-62),  but  are 


multiplied  by  the  normal  reflection  coefficient  T  = 


i  -VH 

l+^^£w 


for  the  wall  with  dielectric 


constant  E1V ,  have  a  different  diffraction  coefficient,  and  phase  factor  due  to  the 
additional  propagation  distance  in  the  vicinity  of  the  comer.  Thus 


Eydif2{x,  z)~  ^ 

“n re~j4K"deJ*  Xl-1)-?-  D^>mr, (])/)cos(P#irv>"2X^ 

-d  [m=odd 

-  S(-OTD((D,/,0/)sin(PHrv)e-2K™4 

m=even  J 

and 


(3-63) 
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(3-64) 


£/'/3M  = 

nJe-j^jenA\  HXy*-= 

2d  {m=»dd 

~  S  (- 1)1  .  6„  )sin(P#„.v>”yK“:  1 


The  received  powers  in  an  averaged  over  the  cross  section  for  vertically  polarized 
sources  in  a  straight  tunnel  and  in  an  L-Bend  tunnel  are  shown  in  Figure  3-16  ~  18.  The 
width  of  both  tunnels  is  4  m  and  the  walls  have  the  dielectric  constant  £,.  =  6  and 
conductance  a  =  0.01 .  Figure  3-16  shows  the  results  for  a  900MHz  source  located  500  m 
from  the  L  bend,  and  centered  in  tunnel.  The  UTD  diffraction  coefficients  [12]  are  used 
for  this  simulation.  The  coupling  loss  of  the  4  m  wide  L-Bend  tunnel  is  about  50  dB  for 
this  case.  Figure  3-17  and  3-18  show  the  received  power  for  300  MHz  and  100  MHz 
sources  located  at  50  m  from  L  bend.  The  turning  loss  for  these  two  frequencies  is  about 
35  dB,  and  30  dB  respectively. 
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Received  Power  (dB) 


Received  Power  in  4m  wide  Tunnel  for  100MHz,  Vertical  Polarization 


Figure  3-16  Received  Power  in  4  m  wide  2-D  Straight  and  L  bend  Tunnel 


900  MHz,  Vertical  Polarization  =  6, a  =  0.01 
The  corner  of  L-Bend  Tunnel  is  located  at  500  m. 


Received  Power  in  4m  wide  Tunnel  for  300MHz,  Vertical  Polarization 


Figure  3-17  Received  Power  in  4  m  wide  2-D  Straight  and  L  bend  Tunnel 
300  MHz,  Vertical  Polarization  ,er  =  6, a  =  0.01 
The  corner  of  L-Bend  Tunnel  is  located  at  50  m. 
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Received  Power  in  4m  wide  Tunnel  for  100MHz,  Vertical  Polarization 


Figure  3-18  Received  Power  in  4  m  wide  2-D  Straight  and  L  bend  Tunnel 
100  MHz,  Vertical  Polarization  ,er  =  6,  a  =  0.01 
The  comer  of  L-Bend  Tunnel  is  located  at  50  m. 
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4.  Propagation  in  Rectangular  3-D  Tunnels 

4. 1  Fields  in  a  Rectangular  3-D  Tunnel  for  Point  Source  Excitation 

The  modal  fields  in  straight  3-D  tunnels  of  rectangular  cross  section  are 
approximately  the  product  of  the  modal  variations  for  2-D  tunnels  [2].  For  the  vertical 
polarization  of  the  electric  field  in  a  3-D  tunnel,  the  power  loss  of  side  walls  is  due  to  TE 
attenuation  and  the  power  loss  of  ceiling  and  floor  is  due  to  TM  attenuation.  Because  the 
spatial  dependence  in  both  transverse  x  and  y  direction  can  be  either  a  cosine  or  a  sine 
variation,  the  modal  field  in  3-D  tunnel  have  four  possible  combination  of  x  and  y 
dependence.  These  combinations  are  :  1)  sine  variation  in  both  x  and  y;  2)  sine  variation 
in  x  and  cosine  in  y;  3)  cosine  variation  in  x  and  sine  in  y;  4)  cosine  variation  both  in  x 
and  y.  To  simplify  the  field  expressions,  we  place  the  transmitter  at  center  of  the  tunnel 
(x'  -  0,  y'  =  0),  so  that  only  the  symmetric  modes  are  excited.  For  vertical  polarization  of 
the  electric  field,  the  symmetric  modes  have  dependence 

Ey  =  £  X  A,,,,  cos(/?„,re.v)cos(/?„7'u  yyjKi"z  (4-1) 

m-odd  n~odd 

For  a  tunnel  of  width  2c/ 1  and  height  2 cU_ 


(4-4) 


9 

The  attenuation  constant  of  the  modes  is  dominated  by  the  term  n~Ew,  which  causes  the 
attenuation  to  be  like  that  of  a  TM  polarization  in  a  2-D  tunnel. 
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The  spatial  dependence  of  the  current  distribution  for  a  point  source  located  at 
y  =  0  can  be  expressed  in  terms  of  the  Fourier  sum 


I0hb(y)=  I0h  +  ^-  X  cos| 

^2  n-even 


f  mi  ^ 


2d , 


(4-5) 


The  cosine  terms  in  (4-1)  can  be  thought  of  as  resulting  from  two  line  sources  with  phase 


,  .  nK 

±7— >' 


variation  e  2d  .  The  phased  line  sources  in  a  2-D  tunnel  produces  total  field  having  y 


dependence  cos 


3 


mi 

\2 d  j 


which  is  close  to  that  of  the  modal  field  in  (2-53).  As  a  result,  the 


modal  field  for  vertical  polarization  in  a  3-D  tunnel  is  given  by 


Ev(.v,y,z): 


r  zIA ' 
^d\d2j 


m~odd  n-odd 


^COPo  A 
V  K'"’"  J 


cos(p,/£*)cos(p/My)e’  JK>" 


(4-6) 


Figure  4-1  shows  the  received  power  in  3-D  tunnel  (4m  wide,  3m  high)  for 
vertically  polarized  900MHz  signal.  The  received  power  is  averaged  over  the  section  at 
1.5  m  high  from  the  floor  and  the  transmitter  is  centered  and  at  the  mouth  of  tunnel.  The 
reference  position  for  normalization  is  10  m  from  the  source  and  it  is  -55dB  at 
reference  point.  The  free  space  received  power  at  this  point  is  -51.5  dB.  The  normalized 
path  gain  for  an  isotropic  point  source  in  a  3-D  straight  tunnel  is  shown  Figure  4-2.  For 
this  calculation,  the  dimension  of  tunnel  is  assumed  to  be  2 d\  =  4  m  and  ld2  =  3  m.  The 
transmitter  is  oriented  vertically  at  the  center  of  the  tunnel  width  and  height.  In  order  to 

define  the  normalized  path  gain,  £  I  is  averaged  over  the  width  of  tunnel  at  1.5  m  high 

from  floor.  This  value  is  then  normalized  to  the  value  of  the  average  at  z  =  10  m.  Because 
the  tunnel  height  is  smaller  than  its  width,  and  because  the  attenuation  due  to  reflection  of 
TM  polarized  waves  at  the  ceiling  and  floor  dominates  over  that  due  to  TE  reflection  at 
the  walls,  the  attenuation  in  the  tunnel  is  greater  than  either  the  TE  or  TM  case  in  a  2-D 
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(gp)- 


tunnel  of  width  4  m.  As  in  the  case  of  a  2-D  tunnel,  lower  frequency  signals  undergo 
much  higher  attenuation. 


Receied  Power  in  4m  X3m  3-D  tunnel  (  900MHz,  Vertical  Polarization  ) 


0  100  200  300  400  500  600  700  800  900  1000 

Distance  (m) 


Figure  4-1  Received  Power  in  3-D  straight  tunnel  (  4m  wide,  3m  high  ) 
for  900MHz,  Vertical  Polarization  (fM.=6,  <7  =  0.01) 
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Normalized  Path  Gain  of  4m  x  3m  Tunnel  for  Vertical  Pol. 


Figure  4-2  Normalized  Path  Gain  in  3-D  straight  tunnel  ( 4m  wide,  3m  high  ) 

for  various  frequencies 
Vertical  Polarization  (e  =6,  <7  =  0.01) 


4.2  Attenuation  of  modal  field  in  3-D  tunnels 

The  attenuation  constant  in  a  3-D  tunnel  that  is  2(7,  m  wide  and  ld2  m  high  is 
defined  by  (4-7)  for  vertical  polarization  and  (4-8)  for  horizontal  polarization  from  (4-6). 


where  eu,  is  dielectric  constant  of  tunnel  walls.  In  Figure  4-3,  we  showed  the  attenuation 
constant  per  100  m  for  dominant  mode  (m=l,n=l)  vs  the  ratio  of  tunnel  width  and  height 
in  3-D  tunnels.  The  tunnel  height  is  assumed  to  be  2 d2  =4  m  high  and  the  operating 
frequency  is  900  MHz  for  this  plot.  The  dielectric  constant  of  tunnel  walls  and  both 
ceiling  and  floor  is  used  eu,=6.  If  the  width  of  tunnel  is  greater  than  the  height,  then  the 

attenuation  of  horizontal  polarization  is  lower  than  vertical  polarization.  Otherwise,  the 
vertical  polarization  gives  lower  attenuation  in  3-D  tunnel.  By  (4-7)  and  (4-8),  the 
polarization  for  low  attenuation  is  determined  by  comparing  the  two  dimension  terms 
with  dielectric  constant. 


The  attenuation  of  dominant  mode  (1,1)  versus  — —  for  3-D  tunnels  having  square 

a 

cross  section  is  shown  in  Figure  4-4.  As  we  discussed  in  connection  with  the  attenuation 
in  2-D  tunnels,  the  attenuation  increases  rapidly  when  the  wavelength  of  given  frequency 
is  close  to  the  tunnel  dimension.  Higher  modes  attenuate  even  more  rapidly.  In  order  to 

quantify  this,  let  A L  =  — - —  be  the  distance  at  which  the  next  higher  mode  (2,1) 

a,i  CX21 

attenuates  by  10  dB  more  than  the  dominant  (1,1)  mode.  The  variation  of  A L  with  the 


parameter 


X2 

d3 


in  a  3-D  tunnel  with  square  cross  section  is  illustrated  in  Figure  4-5.  For 
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A  _? 

example,  in  a  4mx4m  tunnel  and  at  a  frequency  of  900MHz,  — -  -1.42x10  From 

d5 

Figure  4-4  it  is  seen  that  in  just  over  100m  the  (2,1)  mode  is  reduced  by  lOdB  as 
compared  to  the  (1,1)  mode.  At  300MHz,  this  distance  is  only  about  15m. 


Attenuation/1 00m  vs  Width/Height  in  3-D  Tunnels  for  (1,1)  modes 


Figure  4-3  Attenuation  constant  for  dominant  mode  vs  the  ratio  of  the  width  and  height 
in  3-D  tunnels  for  er  =6  (900MHz,  4m  tunnel  height) 
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4.3  Coupling  loss  at  L,  T,  Cross  junctions 


The  coupling  loss  at  junctions  in  3-D  tunnels  will  occur  via  diffraction  at  the 
comers,  as  previously  discussed  for  2-D  tunnels.  If  the  height  of  the  main  tunnel  and  side 
tunnel  are  the  same,  then  the  modes  excited  in  the  cross  tunnel  will  have  the  same  field 
variation  in  the  vertical  direction  as  the  incident  mode.  For  the  lower  modes  the 
diffraction  coefficient  at  the  edge  will  be  nearly  the  same  as  in  the  2-D  tunnel  were  there 
is  no  variation  along  the  edge.  Since  the  coupling  is  nearly  independent  of  the  mode 
number  in  the  vertical  direction,  the  coupling  loss  for  3-D  tunnels  is  the  same  as  for  2-D 
tunnels,  and  is  independent  of  tunnel  height. 

The  coupling  loss  vs  various  tunnel  width  for  L-bend,  T  junction  and  cross 
junction  structure  for  900MHz  vertical  polarization  is  shown  in  Figure  4-6.  The  coupling 
loss  of  cross  tunnels  is  about  3  dB  lower  than  T-junction  tunnels.  The  diffracted  fields  in 
cross  tunnels  have  no  wall  reflection  factor  for  comparing  with  the  T-junction  so  that  the 
coupling  loss  of  T-junction  tunnels  is  higher  than  cross  tunnels.  (See  Figure  3-7  and  3-11) 
The  coupling  loss  of  L-bend  tunnels  is  about  8  dB  higher  than  T-junction  tunnels.  We 
find  the  diffracted  fields  in  T-junction  have  one  more  upper  edge  diffraction  fields  by 
comparing  with  Figure  3-8  and  3-10.  The  incident  angle  for  upper  edge  is  smaller  than 
lower  edge  for  given  same  diffracted  angle  so  that  the  diffraction  coefficient  of  upper 
edge  is  larger  than  lower  edge.  Therefore  the  coupling  loss  of  T-junction  is  smaller  than 
the  L-bend  tunnels. 

If  we  take  into  account  only  the  (1,1)  mode  in  the  main  and  side  tunnels,  then  the 
received  powers  in  the  cross,  T  and  L  junctions  at  y=0  are  respectively 


E,{x,y=0,z)cross  = 


.n 


uine 


4d{d2 


Dr\  +  Dr 


-j  2r/,K, 


- (dC3  +  DCAe  j2<l'K>  1  )e  ■/K,l-f/‘  )cos(p, .v)e" 


AY,-" 


(4-7) 
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£v(x,y  =  0,z)7.  . 


uine 


4  d,d 


junction 

DTi+DT4  —  \Dr'i  +  Dr’iX'e 


j  *  f  ,  r~\  ^  jK |j4rf| 


'T2  _r  ^73 1 


-A'u4^  L- 
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where  uin  is  given  by 


=  £v  (x  =  ±d{ ,  y  =  0,  z  =  L)  ■■ 


- 1  oh  ) 
4  cUd, 


Ku 


cos((3,re(±t/I ))e 


7‘K,.^ 


(4-10) 


and  D,  is  the  UTD  diffraction  coefficient  for  (1,1)  diffracted  mode  from  (1,1)  incident 
mode  at  comer  for  the  width  of  2d\  and  the  height  of  2  d2  tunnel. 

The  coupling  loss  for  (1,1)  mode  can  be  defined  by  the  ratio  of  the  fields  in  main 
tunnel  and  side  arm  and  results  from  comer  diffraction.  By  (4-7)  ~  (4-10)  and  (1,1)  mode 
expression  of  (4-6),  the  coupling  losses  for  three  structures  are  showed  in  Figure  4-7. 
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where  UTD  coefficients  for  coupling  losses  are  listed  in  Table  4-1. 
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(Figure  3-7) 
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Table  4-1  Incident  and  diffracted  mode  angles  for  UTD  coefficient  for  coupling  of  (1,1) 

modes  in  the  main  and  cross  tunnels 


Coupling  Loss  (dB) 


Coupling  Loss  vs  Tunnel  width 
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Figure  4-6  Coupling  Loss  vs  Tunnel  width  by  modal  simulation 

for  Vertical  Polarized  900MHz  source 
(  Dielectric  constant  of  walls  =6  ) 


Coupling  Loss  (dB) 


Coupling  Loss  for  Cross  Tunnel  (Dominant  mode) 


Tunnel  width  (m) 


Figure  4-7  Coupling  Loss  vs  Tunnel  width  for  Dominant  (1,1)  modes 
for  Vertical  Polarized  900MHz  source 
(  Dielectric  constant  of  walls  =6  ) 
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Solid  line  :  Dominant  (1,1)  Modes,  o  :  Multi-modes 
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Figure  4-8  Coupling  Loss  vs  Tunnel  width  for  Dominant  (1,1)  modes  and  Multi-modes 
Vertical  Polarized  900MHz  source  (  Dielectric  constant  of  walls  =6  ) 


4.4  Comparison  of  coupling  losses  with  measurements 


The  simulation  results  of  each  coupling  loss  of  three  junctions  are  compared  with 
the  measurements  by  Sakai  et.  al.  [7].  They  made  three  scaled  tunnels  by  Bakelite  slabs 
sandwiched  between  two  metallic  plates  ( er  =  3.9 ,  e,  =  0.38 ),  the  width  of  L-bend,  T- 

junction  and  cross  tunnels  were  6  cm,  8  cm  and  5  cm,  respectively.  Vertically  polarized 
9.45  GHz  source  was  used,  and  they  measured  signal  at  middle  of  the  other  tunnel  arm. 
To  compare  with  these  measurements,  the  width  of  the  tunnels  is  fixed  as  4  m,  and  the 
operating  frequencies  of  each  tunnel  are  118  MHz  for  cross  tunnel,  189  MHz  for  T- 
junction,  and  141  MHz  for  L-bend  tunnel  with  same  dielectric  constant  of  the  walls  for 
modal  simulations.  The  coupling  losses  by  measurements  were  27  dB  for  cross  tunnel,  32 
dB  for  T-junction  and  more  than  28  dB  for  L-bend  tunnel.  Using  modal  simulations,  we 
have  20  dB  for  cross  tunnel,  33  dB  for  T-junction  and  30  dB  for  L-bend  tunnel. 

Other  measurements  by  Emslie  et.  al.  was  for  comer  loss  in  high-coal  mine  cross 
tunnels  (14  feet  x  7  feet)  [2],  The  dielectric  constant  of  the  walls  was  said  to  be  8,.  =  10. 
The  measured  coupling  loss  of  40  dB,  and  42  dB  for  horizontally  polarized  fields  at  415 
MHz  and  1  GHz.  The  modal  simulation  for  these  measurement  obtained  a  coupling  loss 
of  37  dB  for  415  MHz  and  40  dB  for  1GHz.  The  results  of  these  two  sets  of 
measurements  indicate  the  validity  of  the  modal  expressions  for  prediction  of  coupling 
loss  at  tunnel  junctions.  The  comparison  is  summarized  in  Table  4-2. 


Modal  Analysis 

Measurements 

Cross  Tunnel  [7] 

20  dB 

27  dB 

T-junction  Tunnel  [7] 

33  dB 

32  dB 

L-bends  Tunnel  [7] 

30  dB 

more  than  28dB 

Cross  Tunnel  (415MHz)  [2] 

37  dB 

40  dB 

Cross  Tunnel  (1GHz)  [2] 

40  dB 

42  dB 

Table  4-2  Comparison  of  coupling  loss  with  measurements  [7], [2] 
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4.  Conclusion 


Radio  propagation  in  long  and  narrow  tunnels  can  be  described  using  modal 
analysis.  In  over  sized  tunnels,  there  exists  many  modes  but  only  a  few  lower  modes  are 
dominant  for  propagating  to  long  distance  from  the  transmitter  because  the  higher 
attenuation  constant  of  the  higher  modes.  The  attenuation  constant  in  tunnels  depends  on 


the  mode  numbers,  tunnel  dimensions  and  frequency  through  the  parameter 


and  m 


where  n,  m  are  mode  numbers,  X  is  the  wavelength  and  2 d{ ,  2 r/2 


are  tunnel  width  and  height.  The  attenuation  increases  with  wavelength  (decreases  with 
frequency)  and  mode  number,  and  decreases  with  tunnel  size.  For  tunnels  that  are  wider 
than  they  are  high  (2dl>2d2 ),  the  modes  with  horizontal  electric  field  have  the  least 
attenuation. 


We  have  evaluated  the  modal  coupling  for  tunnels  with  junctions  using  hybrid 
ray-mode  conversion.  The  coupling  mechanism  is  explained  through  mode  diffraction  at 
the  comers  into  the  side  tunnels.  We  have  found  that  the  coupling  loss  is  50  dB  for  L 
bend  junction,  42  dB  for  T  junction  and  37  dB  for  cross  junction  in  4m  wide  for  900 
MHz  signals  that  have  vertical  polarization.  These  coupling  loss  increases  slightly  with 
increasing  tunnel  width. 
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